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The short-time dynamic evolution of an Ising model embedded in an infinitely ramified fractal 
structure with noninteger Hausdorff dimension was studied using Monte Carlo simulations. Com- 
pletely ordered and disordered spin configurations were used as initial states for the dynamic simu- 
lations. In both cases, the evolution of the physical observables follows a power-law behavior. Based 
on this fact, the complete set of critical exponents characteristic of a second-order phase transition 
was evaluated. Also, the dynamic exponent 6 of the critical initial increase in magnetization, as 
well as the critical temperature, were computed. The exponent exhibits a weak dependence on 
the initial (small) magnetization. On the other hand, the dynamic exponent z shows a systematic 
decrease when the segmentation step is increased, i.e., when the system size becomes larger. Our 
results suggest that the effective noninteger dimension for the second-order phase transition is no- 
ticeably smaller than the Hausdorff dimension. Even when the behavior of the magnetization (in 
the case of the ordered initial state) and the autocorrelation (in the case of the disordered initial 
state) with time are very well fitted by power laws, the precision of our simulations allows us to 
detect the presence of a soft oscillation of the same type in both magnitudes that we attribute to 
the topological details of the generating cell at any scale. 

PACS numbers: 05.50.+q, 64.60.Ht, 75.10.Hk, 02.50.-r 



I. INTRODUCTION. 

In the last years, fractal structures with noninteger Hausdorff dimension (du) have attracted the interest of re- 
searchers because these systems, besides serving to model natural materials such as porous rocks, aerogel, etc. 0,0 j 
also offer the possibility of theoretically exploring systems exhibiting critical behavior close to their lower critical 
dimension, i.e., the larger integer dimension in which the system does not exhibit any phase transition at a finite 
temperature. 

The first studies of phase transitions using fractal structures are those of Gefen and co-workers 0, 0, ■ Based 
on renormalization methods, it has been shown that a second-order phase transition at nonzero temperature occurs 
only if the fractal substrate has an infinite ramification order. Moreover, since translational symmetry is a necessary 
condition to proceed with dimensional perturbation Q , the disagreement between the critical exponents determined 
by current methods ((J 0, EJ El EI an d those obtained by continuation of e— expansions to noninteger dimension 
[l2| can be related to the topological features of the fractal structure. 

It is well known that for systems with translational symmetry, the influence of the underlying structure becomes 
negligible at the critical point, i.e., when the correlation length 

is much larger than the cell spacing, and only the dimensionality, the number of components of the order parameter 
together with its symmetry and the nature of the couplings concur to determine the values of the critical exponents 
and the corresponding universality class. However, in fractal systems, where the translational symmetry is replaced 
by scale invariance, the topological details of the generating cell are present at any scale and such universal behavior 
is said to be weak. The critical exponents and the critical temperature depend not only on d,H, but also on the 
connectivity and lacunarity of the fractal El- A direct quantitative study of topological effects has been recently 
published [14j . 

Most of the previously cited studies are based on the same type of fractal, i.e., the Sierpinski Carpet (SC), which 
has an infinite ramification order. Although the same kind of magnetic interaction (Ising model) has been considered, 
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TABLE I: List of critical temperatures and critical exponents reported in the literature for the Ising model in the SC(3,1) 
fractal. The methods used to obtain the data are: RSGR = Real Space Group Renormalization; MC-FSS = Monte Carlo 
simulations and Finite-Size Scaling; MC-RNG = Monte Carlo Renormalization Group; MC-Slope = Monte Carlo simulations 
and the slope method [T(il[ . The index k indicates the generation of fractal used. The Boundary Conditions (BC) used are 
either Periodic (P) or Free (F). 
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these previous studies yield controversial results. Table [H summarizes the list of published results for the case of the 
two dimensional SC with dn — 1.8927, where the generating cell is built by segmenting a square into nine subsquares 
and removing the central one, so that this fractal is termed SC(3,1). In particular, different values have been obtained 
by different authors for the critical temperature (T c , see tabled, which in most cases is lower than T c of the Ising 
model in 2d. Also, a considerable scattering in the data corresponding to the critical exponents (3/v and 7/1/ can be 
observed. It should also be noticed that Monte Carlo simulations and FSS may either predict dn > d e f / 0, IH ° r 
d H < d eff 0. 

One of the reasons explaining the discrepancies reported in the calculation of T c is most likely related to the location 
of the spins on the fractal (see Table |]J, which could be at the vertices [E IE IE E3 or at the center @>IEIEE3 of the 
squares. Consequently the mean number of nearest neighbors per site is not the same. In addition, for the former 
the number of spins as a function of the cell size does not follow a power law and consequently, dn should not be 
expected to enter in the description of the critical phenomena of such systems. Another reason causing discrepancies 
is related to finite-size effects, since scaling analyses reveal very strong scaling corrections for dimensions smaller than 
d = 2 0, IE US- Monceau et al Q have noticed that discrepancies with standard finite-size scaling (FSS) methods 
could also be due not only to the operation of pure (genuine) finite-size effects, but also to a topological contribution 
to scaling corrections They have also pointed out that the hyperscaling law d e ff — 2(3 /i/ + j/u should remain 
valid for d e ff — dn and the susceptibility should follow the expected power-law behavior, which allows to calculating 
the ratio of exponents y/u and the anomalous dimension exponent 77 in a reliable way |15|. In addition, these authors 
have shown the agreement of their results with those determined by Monte Carlo renormalization group techniques 

m ra 

On the order hand, Pruessner et al |lf| have questioned the validity of FSS studies on fractal structures. In fact, 
they have pointed out that each segmentation step represents a new thermodynamic system and cannot be treated as 
a scaled version of the previous one. In order to avoid this shortcoming, the authors have proposed the Slope Method 
|lfj|. The critical exponents obtained using this method suggest that d e // is smaller than du |lOl |. in agreement 
with the 4-e renormalization group prediction. Furthermore, the obtained critical temperature is higher than those 
reported by other authors by using FSS (see Table [Q). 

Another useful approach to obtain independent estimations for the critical temperature and the critical exponents 
is given by the Short- Time Dynamics (STD) method (for available results see rows 1 — 3 of Table ITI) . By 
performing simulations using both the 4th and 5th segmentation steps, Pruessner et al have reported substantial 
differences between the critical temperature estimated by the STD method (see Table HT|) and the data corresponding 
to the FSS approach (see Table [IJ). Furthermore, the exponents vz, 7, and (3 (see also tabic [HJ strongly depend 
on the segmentation step. In addition, Zheng et al. 11] have determined the exponent 6 of the initial increase in 
magnetization which seems to be slightly greater than the figure accepted for the 2d Ising model. However, this 
determination has to be taken with caution because spins are located at the vertices of the fractal and this approach 
is expected to give an inaccurate estimation of T Cl as has already been discussed. 

In view of the scattering of the available data for the critical temperature and critical exponents, the aim of this 
paper is to study the critical behavior of the Ising Model on the fractal structure SC(3,1) using the Short-Time 
Dynamics approach. In order to achieve this goal, we have used a segmentation step bigger than in a previous STD 
study of this system |lfjj . Furthermore, we have confirmed that the time dependence of the magnetization follows a 
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TABLE II: List of critical temperatures and critical exponents reported in the literature for the Ising model in the SC(3,1) 
fractal and obtained using Monte Carlo simulations and the Short-Time Dynamics analysis (MC-STD). PW = Present Work. 
The Boundary Conditions used are Periodic and the index k indicates the generation of the fractal used. 

Ref. T c vz 7 [3 9 z k Spins location 

[11] 2.033(4) - - - 0.211(3) 2.38(4) 7 Vertices 

[10] 1.5266(11) 3.06(11) 1.959(32) 0.1154(29) - - 4 Center 

[10] 1.5081(12) 3.21(15) 2.048(49) 0.120(55) - - 5 Center 

PW 1.4945(50) 3.546(12) 2.22(1) 0.121(5) 0.1815(6) 2.55(1) 6 Center 

Ising 2.269 2.165124] 1.75 0.125 0.19L24] 2.165f24] - 



power-law behavior for times two orders of magnitude larger than in a previous STD study 10]. Additionally, we have 
determined T c and the complete set of critical exponents starting the STD studies with two different initial conditions: 
i) a fully ordered initial state (ground state configuration corresponding to T = 0), and ii) a fully disordered initial 
state corresponding to T = oo. Self-consistency of the results was also carefully checked. On the other hand, we have 
also obtained the value of the exponent 9 related to the initial increase in magnetization at an early time, which have 
not been previously determined for spins placed at the center of the occupied subsquares of the fractal. Finally, the 
dynamic exponent z was evaluated, for first time for this fractal, by means of two independent methods using the 
Binder's cumulant and the critical scaling of time correlation functions. 

The STD approach has shown to be a powerful tool in the study of critical phenomena because the critical exponents 
can be determined before the critical slowing down of the dynamics takes place and they are free from finite size 
corrections, provided that the correlation length £(t) is always smaller than the system size L (£(i) <C L) [l8l Il9|. 
However, in the case that we are interested in, namely the SC(3,1), an intrinsic kind of finite-size effects due to the 
segmentation step k of the fractal can be observed, as it is discussed in details below. Furthermore, the STD analysis 
is particularly useful since the application of FSS to obtain T c is questionable for segmentation steps smaller than 6 
10], which are the typical sample sizes that one is forced to use in practical calculations due to computing limitations 

The outline of this paper is the following: In Section II we describe the magnetic model and the underlying fractal 
structure, in Section III we briefly recall the main features of the Short-Time Dynamics method, and in Section IV, 
the results obtained from the dynamic simulations are presented. Finally, our conclusions are stated and discussed in 
Section V. 



II. THE ISING MODEL ON THE SIERPINSKI CARPET. 



The SC(b,c) is obtained as follows: for each segmentation step (k), a square of length L is segmented into b 2 
subsquares and c 2 subsquares are deleted from the center of the initial square; then the segmentation process is 
iterated on the remaining subsquares. Figure ^ shows a sketch of the SC(3,1), which is used in the present work, 
corresponding to the k = 3 segmentation step. In the limit k — > oo the mathematical fractal SC(b,c) is obtained, and 

the Hausdorff dimension is given by cLh = ~^Tnir • ^ n ^ ne case °f the SC(3,1) the deviation of the mean number of 
near nearest neighbors from that corresponding to the thermodynamic limit (mathematical fractal), using periodic 
boundary conditions, as determined by the transfer-matrix method, becomes negligible for k > 60. 

As mentioned above, this fractal has an infinite ramification order, which implies that the Ising model should exhibit 
a second-order phase transition at finite temperature. Spins were placed at the center of the occupied subsquares. 
Consequently the number of spins increases as a power law of the lattice size, and the exponent is given by <1h- 

The Hamiltonian of the system is given by 



H = -Jj2 a i s i (!) 

where s$ assumes the values ±1, the sum runs over all interacting nearest-neighbor pairs of spins, and the exchange 
coupling constant J is positive (ferromagnetic interactions). 



FIG. 1: Sketch of the Sierpinski carpet SC(3,1) iterated up to k = 3 segmentation step. Spins are placed at the center of the 
filled squares. 



III. SHORT-TIME DYNAMICS APPROACH FOR CRITICAL PHENOMENA. 



According to field theoretical calculations |l7j| , if a magnetic system at high temperature, with a small magnetization 
mo, is suddenly quenched to the critical temperature, it may exhibit a universal dynamic evolution, which sets right 
after a time scale t m i C . It is expected that t m i C should be large in the microscopic sense, but still very small in the 
macroscopic sense necessary for equilibration. This STD approach is free of the critical slowing down since the spatial 
correlation length is still small within the short-time re gime , eyen at (or near) the critical point |17| . 

The fc-th moments of the magnetization are given by [let [19| 



where (3 and v are the order parameter and the correlation length critical exponents, z is the dynamic critical exponent, 
t is the reduced temperature, L is the system size, xq is the scaling dimension of the initial magnetization, and b is a 
scaling factor. For large L, at the critical point r = 0, and for mn <C 1, from the scaling form given by equation (0) 
one derives the initial increase in magnetization, obtaining 18, 19] 



where (...) denotes the averages taken over spin configurations and [...] corresponds to averages taken over different 

samples with equivalent initial conditions. Here, 9 — (^xq — / z, and the scaling function behaves as F(x) ~ 1 for 

x — > and F(x) ~ - for x — ► oo. It should be noticed that the time scale for this initial increase is of the order 

of to ~ m *° . 9 and xq are the exponents of the initial increase and the scaling dimension of the order parameter. 
Since both exponents are related, one of them can be considered as a new no trivial critical exponent [Tsl Il9| . 
Performing simulations for different values of the initial magnetization and extrapolating the results to mo = 0, the 
exponent 9 can be obtained. Other interesting observables are the second moment of the magnetization (M 2 (t)) and 
the autocorrelation (A(t)), that for mo = and r = should behave according to the following power-law scaling 
relationships 



M k (t,T,L,m ) = b- h £M k (b- z t^T,b- 1 L,b Xo mQ) 



(2) 
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and 



A(t) = 



1 N s 
^5>(*) Si(0) 



= 1 



ott-\ with \ = ^ll-e, (5) 



respectively. 

Another important process that can be measured is the dynamic relaxation from a completely ordered state (with 
Too = 1), which corresponds to a ground state configuration at T — 0, to T c . In this case, the magnetization, the 
logarithmic derivative of the magnetization with respect to r, and the second-order Binder's cumulant |24j should 
behave according to 

M(t)oct _ «, (6) 
V T (t)=d T (lnM(t,T)) | T = oct-, (7) 

and 

C/^^-loct^, (8) 
V ' [M(t)f V ; 

respectively. For T ^ T c , but within the critical region, the power-law behavior is modified by a scaling function, 
which for the magnetization is given by M(t~r). This fact can be used to determine the critical temperature from 
the localization of the optimal power-law behavior. 

Summing up, the STD scaling study of a given system performed by starting from two extreme initial states, i.e. a 
completely ordered one and a completely disordered one, is sufficient to determine both the critical temperature and 
the set of relevant critical exponents in a self-consistent fashion 0, 0, |2(| . 

IV. NUMERICAL RESULTS. 

A. Details on the Simulations. 

Monte Carlo simulations of the Ising model on the SC(3,1) were performed for the segmentation step k = 6 
(L = 729 with 262144 spin sites) using periodic boundary conditions, and starting either from an ordered state or 
from a disordered state with zero or a small initial magnetization. In the latter case the initial magnetization (too) was 
obtained from a disordered configuration (of zero magnetization) by flipping a definite number of spins at randomly 
chosen sites. 

In order to implement the time evolution, the system is updated by using the Metropolis algorithm. The well tested 
|2l| Marsaglia-Zanan pseudo random number generator is used throughout the simulations. The time unit, defined 
as a Monte Carlo time step (MCS), involves the update of a number of spins that corresponds to all spin sites of the 
sample. In this way, during one MCS each spin is updated once, on average. Simulations starting from a disordered 
(ordered) state are carried out up to 2000 MCS (200000 MCS). 

The magnetization, the autocorrelation, and the second moment of the magnetization were averaged over a number 
n s of samples with equivalent initial configurations. In addition, the time evolution of the ordered state was also studied 
for the segmentation steps k = 3, 4, 5 in order to apply an FSS method that allowed us to obtain an estimation of the 
dynamic exponent z. In order to estimate the error bars of the evaluated exponents we have used a variant of the 
blocking method (2^ fitting the time dependence of each observable for independent sets of measurements having the 
same statistic. 

B. Simulation of the Dynamic Evolution starting from the ordered state. 



According to our experience, the determination of the critical temperature and exponents is more accurate when the 
simulations start from the ordered state. In fact, in this case the magnetization is large and decreases slowly during 
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TABLE III: List of critical exponents determined from the time dependence of the magnetization (second column, see equation 
©), the Binder's cumulant (third column, see equation JHJ), the logarithmic derivative of magnetization (fourth column, see 
equation 0), and the susceptibility (fifth column). Data obtained starting the simulations from an ordered initial state and 
for the k= 6 generation of the fractal. Slightly different exponents are obtained performing the fits after disregarding different 
initial time intervals (t m in), as listed in the first column. 



tmin(MCS) — till. J_ SL f] j 

" Ltlt \ / y_z Z U_Z V_Z [ 

20 0.03406(6) 0.697(3) 0.285(2) 0.630(3) 0.119(9) 2.21(2) 
100 0.03412(7) 0.693(2) 0.282(1) 0.626(2) 0.121(5) 2.22(1) 
150 0.03413(7) 0.694(2) 0.282(2) 0.627(2) 0.121(9) 2.22(2) 



time evolution and therefore statistical fluctuations are less prominent. Figure [21 shows the decay of the magnetization 
obtained at different temperatures for k — 6. The critical temperature is determined by finding the smallest standard 
deviation from the power law given by equation ©, which yields T c (k = 6) = 1.4945(50), where the error bar is 
assessed by considering the closest pair of temperatures that present noticeable but small standard deviations. This 
is in good agreement with determinations performed by means of Monte Carlo simulations analyzed by using Slope 
method ^(j- Also, acceptable agreement with the value obtained by means of the STD analysis of Monte Carlo 
data is found (see Table [nj. However, a careful inspection of the data shows a systematic decrease in T c when the 
segmentation step is increased, suggesting that our result could be taken as an upper bound. 




FIG. 2: Log- log plots of magnetization versus time obtained starting from ordered initial conditions (mo = 1). The different 
temperatures used and the number of samples considered are also shown in the figure. The inset shows the log-log plot of M(t) 
versus t that gives the best fit of equation © and is assumed to correspond to the critical temperature T c = 1.4945. 

The exponent determined by fitting equation |JBJ, see the inset of Figure|3 is listed in Table ITTT1 (2nd column). Our 
estimation, given by j3/vz = 0.0341(1), is slightly smaller than the figures obtained by using the exponents listed 
in Table CII which correspond to STD studies for k = 4 and 5, namely /3/vz = 0.0377 and (3/vz = 0.0374 [Tol |. 
respectively. So, our finding is consistent with a systematic decrease in (3/vz that is observed when k is increased. 

The logarithmic derivative of the magnetization with respect to r is evaluated by taking the difference between 



7 



the values of M(t) at two temperatures close to T c . The result of this calculation, obtained taking T\ — 1.4795 and 
T 2 = 1.5095, is shown in Figure|3| From this figure it follows that the power-law behavior expected from equation (J7J 
is obtained after t = 70MCS, and the corresponding exponent is listed in Table ITTT1 It should be noticed that using 
T\ = 1.4895 and T2 — 1.4995 one also obtains the same critical exponent (within error bars) but the data are more 
noisy. From the data shown in Table ITTT| it follows that vz — 3.546, which is significantly larger than the exponent 
reported for k = 5, namely vz — 3.21 |lfl| . 




t (units of MCS) 



FIG. 3: Log-log plot of the logarithmic derivative of magnetization versus time obtained at criticality, starting from ordered 
initial conditions (mo = 1). The full line corresponds to the best fit obtained for t > 100 MCS, according to equation J7J|. 

Figure0]shows the determination of the exponent (see the results in Table ITTTf) from the time dependence of the 
Binder cumulant, according to equation ©. Considering that the susceptibility is given by oc U(t) [(AI(t))] 2 oc 
t"* , the exponent 7/1/2 can also be obtained (see inset in Figure0]and Table UTTjt . Notice that the above relationships 
were obtained assuming that the hyperscaling law d e ff = 2(3/v + j/v holds. 

Simulations started from the ordered stated also allow the self-consistent determination of the order parameter 
critical exponent (3 and the exponent 7 of the susceptibility using the exponents listed in the 2nd and 4th columns, 
and in the 5th and 4th columns of Table ITTT1 respectively. The obtained results are also listed in Table UTTI The trend 
of the data for the exponent 7, namely a systematic increase with k, is consistent with the observations reported by 
Pruessner et al. [T^j, which are listed in Table ITT1 for the sake of comparison. However, (3 appears to be less sensitive to 
the change of the segmentation step. It is worth mentioning that our best estimation of the order parameter critical 
exponent given by (3 = 0.121(5) is very close to the exact value corresponding to the Ising model in d = 2 (J3 — 0.125). 
However, 7 is clearly greater in the case of the fractal substrate. 



C. Simulations of the Dynamic Evolution starting from a disordered state 

Figure [S] shows the initial increase in magnetization, observed for different values of the initial (small) magneti- 
zation (too), and obtained after quenching the system to T c when the simulations started from the disordered state 
corresponding to T — 00. Within the time regime considered (20 — 2000 MCS), the magnetization always increases 
and the data can be fitted to a power law with critical exponent 8 7 as expected from equation Nevertheless, a 
soft curvature of the data can be observed for larger times due to the fact that mo is finite and the power law is 
actually expected to hold in the too — > limit. So, in order to determine the critical exponent we performed a fit of 
the data within the time interval 20-100MCS. As can be observed in the inset of figure the exponents show a weak 
dependence on toq. Then, the exponent 8 was evaluated by a linear extrapolation to too = 0, yielding 8 = 0.1815(6). 
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t (units of MCS) 



FIG. 4: Log-log plot of Binder's cumulant versus time obtained at criticality, starting from ordered initial conditions (mo = 1). 
The inset shows the zero field susceptibility (x) versus time. The full lines correspond to the best fits obtained for t > 100 
MCS. 

TABLE IV: List of critical exponents determined from the dynamic behavior of the second moment of the magnetization (second 
column, see equation 10)), and autocorrelation (third column, see equation @). Data obtained starting the simulations from 
disordered initial states with mo = and for k — 6. The estimations of d e ff/z listed in columns 4th and 5th are obtained by 
using the autocorrelation data (3rd column) and the determined value of 6 = 0.1815, and using the exponents listed in the 2nd 
column in combination wich the exponents listed in the 2nd column of Table ITTT1 respectively. 
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So, according to our results, the exponent 9 for the SC[3,1) fractal appears to be slightly smaller than the accepted 
value for the Ising model in d = 2, given by 9 = 0.191 (221. 

Figure El shows that the evolution from a disordered state, with mo = 0, of the second moment of the magnetization 
has a weak dependence on temperature. This shortcoming hinders an independent estimation of T c based on these 
measurements. However, by using the value of T c obtained by means of simulations started from the ordered state, it 
is possible to evaluate the critical exponent of the second moment according to equation J3J , as shown in the inset of 
figure El The obtained value is listed in Table llVl 

On the other hand, the decay of the autocorrelation function (see figure [7J slightly depends on T, allowing us to 
confirm our estimation, namely T c = 1.4945(50), and the corresponding error bars, already evaluated using simulations 
started from ordered configurations. The exponent A = d e tf/z — 9 evaluated by fitting the data at criticality (see 
inset of figure 0) is also listed in Table ITVl It is worth mentioning that by inserting the value of 9 already determined 
in the exponent of the autocorrelation function, one can also calculate d e ff/z, as listed in the 4th column of Table II VI 
The obtained results are in full agreement with the determination performed by starting simulations from ordered 
states, sec Table Hill 
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FIG. 5: Log-log plots of magnetization versus time obtained at criticality, starting from disordered initial conditions slightly 
modified to obtain different values of the initial magnetization mo. Data corresponding to k = 6 and different values of mo, 
which from top to bottom are: 0.06, 0.05, 0.04, 0.03, 0.02, 0.0f5 and 0.01 , respectively. The full lines correspond to the best fits 
obtained for 20 < t < 100 MCS, according to equation Data obtained by averaging over 4000 — fOOOO different samples, 
depending of mo. The inset shows the dependence of 6 on the initial magnetization mo that allowed us to extrapolate the 
exponent 9(mo — * 0) = O.I8f5(6). More details in the text. 

Also, inserting the exponent (3/vz determined by means of simulations started from the ordered state (see the 
2nd column of Table in the expression of the exponent of the second moment of the magnetization, given by 
d e ff/z — 2/3/vz, one can obtain an additional estimation of d e ff/z, as listed in the 5th column of Table HV1 which 
is slightly smaller than the estimations already performed by means of different procedures. We attribute this small 
difference to the propagation of errors in the evaluation of exponents by combining results from different measurements. 

D. Determination of the dynamic exponent z. 

Although the application of FSS techniques to the evaluation of critical exponents in systems with fractal structure 
is questionable, as we have already mentioned in the introduction, we used it at this point in order to obtain a 
first estimation of dynamic exponent z. However, we will perform a critical analysis of the obtained results and 
subsequently we will perform a second (more accurate) estimation of z. We carried out simulations up to 10 4 MCS 
for segmentation steps k = 3, 4, 5 and 6 (L = 27, 81, 243 and 729, respectively) and determined the dynamic exponent 
z using an FSS analysis of the Binder's cumulant. In fact, right at T c , the dynamic exponent z can be determined 
from the Binder's cumulant according to the following scaling relation 

U(t,L 1 ) = U(tb z ,L 2 ), (9) 

where b = . Figure|Slshows the data collapse obtained when the time scale of the system of size Li size is rescaled by 

a factor (ff) • The results obtained by rescaling lattices of sizes 81/27, 243/81, and 729/243 are z = 2.76(2), 2.65(1) 

and 2.60(3), respectively. So, we observed a systematic decrease in z when the segmentation step of the fractal is 
increased, and consequently it is no longer valid to set a single value of the dynamic exponent for all segmentation 
steps. This observed behavior is similar to an observation reported previously Q where the fixed point intersection 
of the Binder's cumulant for different sizes was replaced by a sequence of intersection points occurring at "effective" 
critical temperatures, while the actual critical temperature was defined as the limit for k — > oo. 
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FIG. 6: Log-log plots of the second moment of magnetization versus time, obtained at different temperatures (listed in the 
figure), starting from disordered initial conditions with mo = 0. Data obtained by taking k = 6. The inset shows the dependence 
of M 2 on time obtained at criticality (T c = 1.4945), where the exponent d " if is obtained by fitting the data according to 
equation Q. See also Table HVI Data obtained by averaging over 8000 different samples. 



Summing up, the FSS method applied to the dynamic behavior of Binder's cumulant only allow us to establish an 
upper bound to the dynamic exponent given by z = 2.60(3). 

In order to obtain an independent estimation of the dynamic exponent z that doesn't involve calculations with 
segmentation steps smaller than 6 we study the scaling behaviour of the time correlation function: 



C(r,t) 



N 



^Y^Si(t) s l+r {t) 



(10) 



where i + r indicates a site displaced by r lattice spacings relative to site i. Our purpose is to study the onset of 
correlations between spins when an initially completely disordered system (T = oo) has been quenched to T — T c . 
Conventional critical scaling implies the following scaling form for C (r, t) : 



C(r,t) = r-^- 2+ ^f c (r/m)- (H) 

Assuming that the hyperscaling relation given by d e ff — 2(3 jv + j/v holds for this system, and using 77 = 2 — j/v, 
we may replace d — 2 + r\ in equation l|ll|) by 2/3 jv. As 2(3 /vz has already been obtained directly in the simulations 
from the decay of the magnetization from the ordered state (equation ©), and £(t) is expected to behave as t x l z , we 
may plot 



M>C{r,t) vs (12) 

and look for the value of z that make the curves to collapse. This procedure has been applied for Humayun and Bray 
|25j to obtain z for the Ising model for d = 2. 

In figure El we show plots of C(r,t) as a function of t obtained for different values of r ranging from 4 to 23. The 
inset shows the best collapse of the curves obtained for z = 2.55. This value was obtained performing a fit of the 
scaled data to a 4-parameter function given by 
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FIG. 7: Log-log plots of autocorrelation versus time, obtained at different temperatures (listed in the figure), starting from 
disordered initial conditions with mo = with k = 6. The inset shows the dependence of autocorrelation on time obtained at 
criticality (T c = 1.4945), where the exponent A = -^-—d is obtained by fitting the data according to equation ©. See also 
Table HV1 Data obtained by averaging over 8000 different samples. 



that we have empirically found to fit the data quiet well. We have also checked that other functional dependences for 
f(x) that also fit the data only modifies the value of z by less than 0.01. The obtained value for k — 6, i.e. z — 2.55(1) 
is consistent with the trend observed from the FSS analysis of Binder's cumulant and sets our upper bound for the 
dynamic exponent. 

The relationship d e ft/z = 0.693 has been determined quite accurately by using three different kinds of measure- 
ments: the Binder's cumulant, the slope of the second moment of the order parameter and the autocorrelation function 
(see Tabic IhT|) . So, taking z — 2.55 the effective dimension becomes d e ff ~ 1.77 , i.e. a figure that is noticeably 
smaller than the Haussdorf dimension dn = 1-89 of the SC(3,1). It should be noted that previous estimations us- 
ing FSS studies gave values of d e ff very close to da (for a compilation of published results see table [!}. However, 
Pruessner et. al jlOI | . have estimated d e ff ~ 1.7, also smaller than the Haussdorf dimension. 

On the other hand, using the value 1/vz — 0.282 determined from the slope of the logarithmic derivative of the 
order parameter (see equation (J7J) and z — 2.55, our estimation for the lower bound of the correlation length exponent 
becomes v ~ 1.39, which is significantly larger than the exact value of the Ising model given by v = 1. 



It is shown that the Short-Time Dynamics approach is a powerful method for the study of the critical behavior of 
the Ising Model embedded in a fractal structure, where the translational symmetry is changed for the scale invariance. 
This method allows us to obtain a self-consistent determination of the critical temperature and the complete set of 
the critical exponents. This self-consistence is achieved by using three different initial conditions for the study of the 
dynamics. 

The critical temperature determined in the present work for k — 6 (an upper bound for this system) is in agreement 
with the value reported by Prussener et al. |10| . which was obtained from equilibrium measurements. We note 
that these authors have used the Slope method, a procedure that is free of the finite size effects involved in FSS 
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FIG. 8: Scaling plots of the second-order Binder's cumulant obtained using data corresponding to adjacent pairs of segmentation 
steps (ki,ki+i) with sizes (Li, Li+i), respectively. The crosses (solid lines) correspond to systems of sizes Li+i (Li). The time 
scale for the latter is rescaled by a factor (Li+i/Li) z in order to achieve superposition to the former; in each case the exponent 
z is taken as a fitting parameter. Data obtained by averaging over 720, 17000, 47000 and 70000 different samples, for k — 6, 5, 4 
and 3, respectively. 

calculations. A critical discussion of the values reported by other authors employing different techniques is presented 
in the discussion of the results listed in Tables HI and ITU 

The exponent 9 of the initial increase in M(t) determined for the segmentation step k — 6 and extrapolated to the 
too ~ y limit is slightly smaller than the value corresponding to the two dimensional Ising model (see Table QJ. Our 
result is also in agreement with that reported by Zheng et al. [r7J. obtained by locating the spins at the vertices of 
the SC(3,1) and for the same segmentation step. In this way, these results suggest that the fractal structure does not 
significantly affect the exponent 9. 

The dynamic exponent z has been obtained by means of two independent measurements. Binder's cumulant method 
allow us to determine the decreasing trend of z when k is increased, so we obtained an upper bound given by z = 2.60. 
Further analysis of correlation functions allow us to improve the estimation of such an upper bound which is given 
by z = 2.55. 

Our value for the exponent (3 = 0.121 of the order parameter for the segmentation step k — 6, is also in agreement 
with the trend of the results reported by Pruessner et al. 01 f° r segmentation steps k — 4 and k — 5 (see Table 
II), suggesting that our estimation can be taken as a lower bound. This value is only slightly smaller than the exact 
exponent corresponding to the 2d Ising magnet, namely (3 = 0.125. So, we should not disregard the possibility that 
for k — > oo the order parameter critical exponent may adopt the same value for both systems. 

On the other hand, our estimations of the exponents 7 = 2.22 of the susceptibility and v = 1.39 of the correlation 
length, are significantly larger than those obtained for the 2d Ising system, namely 7 = 1.75 and v = 1, respectively. 
Observing the trend reported by Pruessner et al. [TJj (see Table II) it is expected that our estimation for 7 could be 
taken as a lower bound. Also, v may be taken as a lower bound because it is evaluated from the measurement of vz 
that involves our estimation of an upper bound for z. 

The relationship d e ff/z = 0.693 has been determined quite accurately, then taking z = 2.55 we conclude that 
d e ff ~ 1.77 < d H . 

Finally we would like to remark that due to the huge statistic achieved in the evaluation of the dynamic properties, 
a soft oscillation around the power-law decay of the magnetization was observed at criticality (see figure [2J . The 
same oscillation was also observed for the disordered initial state in the behavior of the autocorrelation (see figure [TJ. 
The oscillation can clearly be detected by subtracting the fitted power law from the actual data, as shown in figure 
ITUl This oscillation is very nicely reproduced in both measurements (up to t ~ 2 x 10 3 in figure ITU)! and, to the best 
of our knowledge this is the first evidence reported about this interesting behavior of the dynamic properties of the 
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FIG. 9: Time correlation functions C(r,t) for r ranging from r = \/T7 (upper curve) to r = (1 ower curve) obtained for 



the SC(3,1) and k — 6. The inset shows plots of the scaled correlation y — r\"*> C(r,t) as a function of the scaled variable 
x — r/t x ' z , for a = 2.55 and 22 =0.0682. We have taken the value of 22 f r0 m the decay of M{t) obtained in this work and 
adjusted the value of z that gives the best collapse of the curves. For z — 2.50 and z — 2.60 the curves show noticeable 
deviations from the collapsed form (not shown here for the sake of space). 

Ising model in a fractal substrate. We have clear signs that the observed oscillations are related to the topological 
properties of the fractal lattice. A more detailed investigation of these oscillations for the SC(3,1) and other fractals 
will be published elsewhere 26]. 

We would finally like to remark the self-consistency of our results obtained upon the application of the STD method 
to the study of the critical behavior of the Ising model on fractal structures. However, it should be recognized that 
there are still discrepancies in the values of the critical exponents when STD-results are compared to those obtained 
using standard finite-size scaling of equilibrium data. For example, in contrast to our data, recent FSS results of 
Carmona et al and Monceau et al. lead to an effective dimension almost equal to the Hausdorff one. So, we 
conclude that the origin of the discrepancies may be related to the fact that the critical behavior of the Ising magnet 
on fractal substrata is very particular, since it is linked to the dependence of most physical observable upon the 
number of iteration steps of the structure. 
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